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Abstract 

Our study is dedicated to the probabilistic representation and numerical approximation 
of solutions to coupled systems of variational inequalities. The dynamics of each component 
' of the solution is driven by a different linear parabolic operator and suffers a non-linear de- 

■ pendence in all the components of the solution. This dynamics is combined with a global 

structural constraint between all the components of the solution including the practical ex- 
ample of optimal switching problems. In this paper, we interpret the unique viscosity solution 
to this type of coupled systems of variational inequalities as the solution to one-dimensional 
\ constrained BSDEs with jumps introduced recently in [6]. In the spirit of 0], this new repre- 

Q,^ ' sentation allows for the introduction of a natural entirely probabilistic numerical scheme for 

^\ . the resolution of these systems. 

o 
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1 Introduction 

Pardoux and Peng [12] developed the theory of Backward Stochastic Differential Equations, 
providing a probabilistic representation of solution to quasilinear parabolic PDEs. Coupling 
the diffusion process with a pure jump process, Pardoux Pradeilles and Rao [13] extend this 
representation to systems of coupled semilinear PDEs with different linear differential operators 
on each line. Introducing restrictions on the domain of the backward process. El Karoui et al 
[5] cover the class of variational inequalities. Constraining instead the jump part of the solution, 
Kharroubi, Ma, Pham and Zhang [llj allow to consider quasilinear variational inequalities. 

The focus of this note is to extend this type of Feynman-Kac representation to the more 
general class of coupled systems of quasilinear variational inequalities, arising for example in 
optimal impulse or switching problems. We will typically consider systems of PDE of the form 

Cvi - fi{.,{vk)i<k<rn,(yj D:,Vi) A min hij{.,Vi,Vj,aJ D^Vi) = 0, (1.1) 



dt 



l<j<m 



on I X [0, T) X R'^, with terminal condition Vi(T, .) = gi on T x M"* , (1.2) 
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where, for any i G 2 := {1, . . . , m}, is a linear second order local operator 

Cvi{t,x) := bi{x) ■ DxVi{t,x) + ^tT{aiaJ{x)Dlvi{t,x)) , (1.3) 

and 6, a, f, h and g are Lipschitz continuous functions. As observed by [2], this PDE appears in 
the resolution of optimal switching problems as well as stochastic target problems with jumps. 
The major difficulty tgHgs in the coupling between all the components {vi^i<.rn 

of the solution 

and the use of different linear operators at each line. When m is high, the numerical resolution 
of (ll.ip - (ll.2p by classical PDE approximation methods is very tricky and highly computational, 
and we intend to provide here a probabilistic representation to (ll.ip - ()1.2p leading to an efficient 
probabilistic numerical scheme. Whenever b and a are independent of the regime i £ I and the 
constraint functions are of the form hij : {.,yi,yj, .) ^ yi — yj — Cij, Hu and Tang [10] interpret 
the vector solution to (|l.ip - (|1.2p as a multi-dimensional BSDE with terminal condition and 
oblique reflections. The challenging derivation of a convergent numerical approximation for this 
type of BSDE is of great interest and is currently under study. The approach of this paper relies 
instead on a recent reinterpretation of obliquely multi-dimensional reflected BSDEs in terms of 
one-dimensional constrained BSDE with jumps, as introduced in [6]. The idea is to consider, 
as in [13], a random regime driven by a pure jump transmutation process, allowing to retrieve 
simultaneously some information concerning all the components of the solution. 

Being given a (i-dimensional Brownian motion W and an independent Poisson measure fj, on 
M+ X T, we consider, for any initial condition e := {t,i,x) G [0,T] x X x =: the unique 
I X M'^-valued solution of the SDE: 

f Is = i + It - Ir~ )Kdr, dj) t<s<T (14) 

\ Xs = x + f'b{Ir,Xr)dr + f'a{Ir,Xr)-dWr ' "'^ ■ ^'^ 

Formally, given a smooth solution v to (ll.ip - ()1.2p . the process Y := vj£{.,X^) satisfies 

Yt = 9r^{X^) + rfiiiXt, + Us, Zs)ds + Kt - Kt - [\ ■ dW^ - Cf Us{j)^i{ds, dj) (1.5) 
Jt Jt Jt Ji 

on [0,T], where we denote Zs := (tX {X^)DxVj<^_ (s,Xf), Us{.) := vX-s^X^) — vj^_ (s,Xf), and 

■■= /o*[ — -C^^vie^ - fie{., {vk)i<k<m, crJeD^vie){u, X^)]du. Slucc V satisfies ([Til, we 
expect the following constraint to be satisfied: 

hie^{XI,Ys-,Ys- + Usij),Zs) > 0, jel, t<s<T. (1.6) 



The BSDE (jl.Sp combined with constraint (11. 6p enters into the class of constrained BSDEs 
with jumps and admits a unique minimal solution under mild conditions on the coefficients. We 
reinterpret the y-component of the solution as the unique viscosity solution to the coupled system 
of variational inequalities p.ip - (|1.2p . This new Feynman-Kac representation is meaningful to 
the BSDE literature since: 

• It extends the results of [11] to more general constraints and driver functions depending 
on U, allowing a strong coupling between the dynamics of the value function components 
and gives a minimality condition in some particular cases. 
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• It enlarges the conclusions of Peng and Xu [T^ derived in the no-jump case. 

• It offers a PDE representation to reflected BSDEs with interconnected obstacles introduced 
in [9], since they relate directly to constrained BDSE with jumps, see [6]. 

• It generalizes the use of diffusion-transmutation process in [13] to systems of variational 
inequalities. 

This representation leads to a natural probabilistic algorithm for the resolution of (jl.l|) - ()1.2p . 
The constrained BSDE with jumps is replaced by a penalized BSDEs with jumps, which is 
approximated by the dicrete time scheme studied in [3] and [7]. This leads to a convergent 
numerical scheme based on time discretization, Monte Carlo simulations and projections. 

The rest of the paper is organized as follows: In Section 2, we discuss existence, uniqueness, 
penalization and give a minimality condition for constrained BSDE with jumps (|1.5p - (|1.6|) . Sec- 
tion 3 presents the viscosity properties and the last section details the numerical approximation. 

Notations. Throughout this paper, we are given a finite horizon T and a probability space 
(r2,^,P) endowed with a (i-dimensional standard Brownian motion W = (Wt)t>o, and an inde- 
pendent Poisson random measure /U on x X, with intensity measure \{di)dt for some positive 
finite measure A on X := {1, . . . ,m}. We denote E := [0,r] x T x W^. For a smooth function 
if : [0, T] X M"^ X T — > M, Dxf and D'^ip denote resp. the derivative of ip w.r.t. i, the gradient 
and the Hessian matrix of (p w.r.t. x. The dependence in w G is omitted whenever explicit. 



2 Constrained Forward Backward SDEs with jumps 

We present in this section the constrained Forward Backward SDEs with jumps and recall the 
existence and uniqueness results of [6]. We discuss the correspondence between the value function 
associated to Y and the U component of the solution. Under additional regularity of the value 
function, we provide a Skorohod type minimality condition for the considered BSDE. 

2.1 Existence and uniqueness of a minimal solution 

As discussed above, the forward process is a transmutation-diffusion process composed by a pure 
jump process / and a diffusion without jump X whose dynamics depends on /. For any initial 
condition e := {t,i,x) € E, {P,X^) is the unique solution to ()1.4p . starting from (i,x) at time t. 



For any initial condition e G E, a solution to the constrained BSDE with jumps is a quadruplet 
{Y% Z% C/^ K^) G 52 X L2, X L? X A2 satisfying where, 

• is the set of real valued ^-adapted cadlag processes Y on [0, T] such that 



||y||52 := E [supo<^<T ' < 00, 

• is the set of predictable M"^- valued processes Z s. t. ||Z||l2^ := E 

• L? is the set of "P cr[I) measurable maps U : ^ x [0,T] x I - 



\U\ 



Q I 1 2 dv 



such that 



E 



< 00, 



< 00, 



is the closed subset of composed by nondecreasing processes K with Kq = 0. 



3 



Furthermore, {Y, Z, U,K) is referred to as the minimal solution to (jl.5p - (jl.6p whenever, for any 
other solution (Y', Z' , U' , K'), we have Y <Y' a.s.. In order to ensure existence and uniqueness 
of a minimal solution to (jl.5p - (jl.6p for any initial condition, we impose the following assumptions. 



(HO) The following holds: 

(i) For any G T^, fi, gi and /ijj- are Lipschitz functions with linear growth. 

(ii) The function hij{x, y, ., z) is non-increasing for all (i, x, y,z,j) G X x M"^ x M x M'^ x I. 

(iii) There exist two constants Ci > C2 > —1 and a measurable map 7 : J X M*^ X M X M'^ 
[M^]2 X J ^ [C2, Ci] such that, for any (i, x, y, z,u,u') e I x M.'^ x M. x x [M^]2, 



fi{x,y + u,z) - fi{x,y + u,z') < J^{uj - u'j)j{i,x,y, z,u,u' ,j)X{dj) . 

(HI) For any e = {t,i,x) G E, there exists a quadruple {Y'^ , Z'^ , If^ , K'^) G 5^ x x L? x 
solution to (|1.5p - (|1.6p . with Yf = vi^{t,Xf), for some deterministic function v satisfying 
\vi{t,x)\ < C(l + \x\) on E. 

We provide in Remark 13.21 a more tractable sufficient condition under which (HI) holds. 

Theorem 2.1. Suppose (HO)-(Hl) holds. For any e := {t,i,x) G E, there exists a unique 
quadruple (Y^ , Z'^ ,U'^ , K^) £ x LL x L? x minimal solution to (|1.5p - (jl.6p . and 



Vi 



{t, x) 1-^ Y^'^'^ defines a deterministic map from E into 



Proof. This result is a direct application of Theorem 2.1 in [6]. □ 
2.2 Related penalized BSDE 

For any initial condition e £ E and n G N, we denote by (ye,n^ ^e,n^ jye,n-j solution to the 
following penalized BSDE with jump 

Yt = gie{X^) + j^"^ fie{XI,Ys + Us,Zs)ds- I^Us{j)fi{ds,dj)- j^^ Zs-dW 

+ n j^[hi._,,{Xt,Ys-,Ys-+Us{j),Zs)]-\{dj)ds, 0<t<T. (2.1) 

Under (HO), following the arguments of [T], we verify that there exists a unique solution to (|2.ip 
and introduce K^'"- := fj[hje^ j[X^ ,Y^1^ jY^l"" + Us'^{j), Zp")]" X{dj)ds , for any e £ E and 
n G N. Furthermore, it converges to the solution of (|1.5p - (|1.6p . 



Proposition 2.1. //(HO)-(Hl) holds, (y^'*") 

„gN converges increasingly to Y^, for any e £ E. 
Additionally, If the process Y^ is quasi-left continuous in time, we have 

||ye_ye,n|| ^ ^||^e_^e,n|| ^ ^ ||fje _ fje,n|| ^ ^ ||^e _ ^e,n|| ^ ^2.2) 

Proof. Fix e £ E and observe from Proposition 2.1 in [6] that y^'" converges increasingly to 
y^. Since // is a Poisson measure, the process y^'" is quasi-left continuous. If Y^ has the same 
regularity, the predictable projections of Y^ and y^'" are simply given by (5^1 )t and {Y^_l^)t 
leading to YfL = lim„_^oo We deduce from the weak version of Dini's theorem, see ^ p. 



4 



202, that y^'"- converges uniformly to Y'^ on [0, T], and the dominated convergence theorem leads 
to ||y*^ — y*^'"!! — > O. Combined with standard estimates of the form 



\Z' 



e,n+'p 



U' 



e,n||2 



_|_ II ^e,n+p _ j^e,n II 2 ^ (^||ye,n+p _ y 



e,n||2 
S2^ 



this implies that the sequences (^"), (C^") and (i^") are Cauchy and hence convergent. 



□ 



Remark 2.1. Under the additional Assumption (H2) below, (fj)igi interprets as the unique 
viscosity solution to p.l|) - (|1.2|) . see Theorem 13.21 In this case, (ui)iex is continuous, Yt = 
vi^{t,Xt) is quasi left continuous and Proposition 12.11 holds. 

We denote by {v^)nen the sequence of deterministic functions defined by : e £ E y/'"" 
and we shall use indifferently the notation v"{t,i,x) or v'!^{t,x), for {t,i,x) £ E. Under (HO)- 
(Hl), we know from Proposition 12. II that v is the pointwise limit of (f")neN- 



2.3 Representation of U and minimality condition 

Proposition 2.2. Let (HO)-(Hl) hold. For any e £ E and stopping time valued in [t,T], we 
have Yq = vje(^9^ X^) , and the process U represents as 



u!{j) 



Vj{s,Xl)-vje_{s,Xl). 



j £l, t<s<T 



(2.3) 



Proof. According to Proposition 12.11 we simply need to provide similar representations for 
the penalized BSDE (j2.ip . Fix e £ E. For any stopping time 6 valued in [t,T], uniqueness of 
solution to (j2.ip and the Markov property of {P, X"^) directly lead to Yg'^ = {d, Xg). Denoting 
i7|'"(j) := v'^is, XI) - vf, (s, X|), for j £l and < s < T, we deduce from ([Ml) that 



Up^{j)fi{ds,dj) 



Yf 



Yf 



U^'-U)f^ids,dj) 



< s < T. 



Therefore E 



lo — Up"' (j))"^ \{dj)ds = and the proof is complete 



□ 



Under an extra regularity assumption on the function v satisfied under Assumption (H2) 
below, the previous representation leads to a Skorohod type minimality condition for (jl.5p - (|1.6p . 

Corollary 2.1. Let (HO)-(Hl) hold. Suppose {vi)i<^x is continuous and the function h does not 
depend on z. Then, for any e £ E, the minimal solution (Y^, Z^, U^, K^) satisfies 



T 

min 

t ie^ L 



/i,e_,,.(x|,y;_,y/_ + c/|(j)) 



dKl 



0. 



(2.4) 



Proof. Fix e £ E. Since {vi)i<zx is continuous, the process Y'^ inherits the quasi-left continuity 



of (/^X'^). Combining and Proposition O leads to maxjgi - f/^'"(i)||^: 

deduce from (j2.2p and Lemma 5.8 in [8], which also holds for caglad functions, that 



mm 

t JG^ 



hi^_,^{XtXl\Y^l- + Ur{j)) 



dKt 



T 

min 



hi^_,{xt,Y:^,Y:^+ut{j)) 



O.We 



die 



Since ^'^ miujgi 



t iGX L 

hie_j{XI, Y^l", Y^f + ;7s'"(j)) dKT < and ^ holds, we get □ 
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3 Link with coupled systems of variational inequalities 



In this section, we interpret the minimal solution to (|1.5p - (|1.6p as the unique viscosity solution 
to the PDE (jl.ip - (jl.2p . thus generalizing the representation derived in [llj . |13j and |14j . 

3.1 Viscosity properties of the penalized BSDE 

The penalized parabolic integral partial differential equation (IPDE) associated to (j2.ip is nat- 
urally defined, for each n G N, by 

— 0^ - -CVi - fii-, Dx^i) - n [hij{.,ipi,ipj,aj D^ipi)]' Xidj) =0 

on [0, T) xR'^ X I, and Vi{T, x) = gi on I x R'^ , 

where C is the m-dimensional Dynkin operator associated to X and defined in (jl.3p . Since the 
penalized BSDE enters into the class of BSDE with jumps studied by Pardoux, Pradeilles and 
Rao [13], we deduce the following Feynman-Kac representation result. 

Proposition 3.1. Under (HO)-(Hl), the functions (v")n are continuous viscosity solutions 
to ([XT]) . Indeed, for any n G N, u"(r, .) = g and, for any ii,t,x) G J x [0, T) x M"^ and 
(f G C'^'^([0, T] X M*^) such that {t,x) is a null global minimum (resp. max.imum) of {vf — if), 
we have 



X 



{t, x) > {resp. < ) 0. 



Proof. Fix n G N. The continuity of follows from similar arguments as in the proof of 
Lemma 2.1 in [13]. According to the representations detailed in the proof of Proposition 12.21 
the viscosity property of v"' fits in the framework of Theorem 4.1 in [13], up to the comparison 
theorem for BSDE, which is replaced by Theorem 2.5 in Royer [15]. □ 



3.2 Viscosity properties of the constrained BSDE with jumps 

Formally, passing to the limit in (j3.ip when n goes to infinity, we expect v to be solution of (jl.ip 
on [0,T) X M*^ X I. As for the boundary condition, we can not expect to have v{T~, .) = g, and 
we shall consider the relaxed boundary condition given by 

min 



Vi - gi , mmhij (.,Vi,Vj,a D^Vi) (T ,x) = on T x M . (3.2) 
jGX \ / J 

Remark 3.1. In the particular case where the driver function / is independent of {y,z,u) and 
the constraint function is given by hij : {x, y,y + v, z) i— > —Cij — v with c a given cost function, 
we retrieve the system of variational inequalities associated to switching problems 

- Cvi - fi, min [vi - Vj - dj] 1 = , on [0, T) X M'^ X J , (3.3) 

Ot j£X J 



mm 



mm 



g^, miujej [vi - Vj - aj]] (T", .) = , on R'^xl . (3.4) 



Thus, if (j3.4p satisfies a comparison theorem, v{T~ , .) interprets as the smallest function grower 
to g satisfying (j3.4p . In particular, we retrieve the terminal condition v{T~ , ■) = g proposed by 
[lOj when the terminal condition g satisfies the cost constraint. 
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In order to define viscosity solution of (jl.ip - (j3.2p . we introduce, for any locally bounded 
vector function (nj)jgx on [0, T] x M'^ its lower semicontinuous and upper semicontinuous (Isc 
and use in short) envelopes and u* defined, for {t, x) G [0, T] x M'^, by 

?x*(t,x)= liminf u{t',x'), and u*{t,x) = limsup u{t',x'). 

{t',x')^{t,x),t'<T {t',x')^{t,x),t'<T 

Definition 3.1. A vector function {ui)i(zx, Isc (resp. use) on [0,T) x W^, is called a viscosity 
supersolution (resp. subsolution) to (jl.ip - ()3.2p if, for each {i,t,x) G X x [0, T] x and if G 
C"^'^([0, T] X M'^) such that {t,x) is a null global minimum (resp. maximum) of{ui — ip), we have, 

dip 

ift<T, min - ^ - C'(p - fi{., {uj)j(zx,aj D^ip),mmhij{.,Ui,Uj,aJ D^tp) {t,x) > {resp. <)0, 

Ot j£X 



if t = T, min 



Ui-Qi, m.inhij{.,Ui,Uj,ai D^(p) (T, x) > ( resp. <) 0. 
jei J 



A locally bounded vector function {ui)i^x on [0, T) xM*^ is called a viscosity solution to (|l.ip - ()3.2p 
if u-t and u* are respectively viscosity supersolution and subsolution to (jl.ip - (l3.2p . 

Theorem 3.1. Under (HO)-(Hl), the function v is a (discontinuous) viscosity solution to 

(inD-dssi). 

Proof. First, following the lines of the proof of Lemma 3.3 and Remark 3.2 in [TT], standard 
estimates on the penalized BSDE (12. 1[) lead to 



E 


sup 






Lte[o,T] J 





\gj.(X;r)\^+ r \Xt\^ds+ sup \vi.{s,Xt)\'' ], eGE. 

Jt sG[0,T] 



Combining Fatou's lemma with standard estimates on X and linear growth conditions on g and 
V, see (HI), we get that sup^g^.T] ^ ^'(l + lxp) with C > 0. Thus, v is locally bounded. 

We observe that the viscosity property of v in the interior of the domain is based on the same 
arguments as the one presented in Theorem 4.1 of [11]. The only difference relies on the more 
general form of the coefficients / and h, which is not a relevant issue since they are continuous. 
In order to alleviate the presentation of the paper, we choose to omit it here and detail only the 
viscosity property (j3.2p on the maturity boundary. 

(i) Let first consider the supersolution property of to (|3.2p . Let (i,xo) G X x M'^ and G 
C"'^'^([0,r] X W^) such that (T, xq) is a null global minimum of — if). Passing to the limit 
the viscosity properties of the penalized BSDE, we get 

vainhij{x,[v^(\i,[v^]j,a'^ Dxip)]{T,XQ) > 0. 

Furthermore v^{T, .) = g, n £ N, so that the monotonic property of the sequence of continuous 
functions {v^'')nen leads to v^{T, .) > g. Therefore v=k is a viscosity supersolution to (13. 2p . 

(ii) We now turn to the subsolution property of v* . Let reason by contradiction and suppose the 
existence of {i,xo) G X x M"^ and G C^''^{[0,T] x M°') such that 

= {vt - 'p){T,xo) = max (v*- If), (3.5) 
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and mill ip — gi , mmj^2hij{.,{p,Vj,aJ Dx^p) {T,xq) =: 2e > 0. The regularity of v* , ip and 
as well as the monotonic property of h lead to the existence of an open neighborhood O of 
{T,xq) G [0,r] X W^, and T,r > such that for aU {t,x,r],r]') G O x (-T,T) x -8(0, r), we get 



mm 



if-r]- gi , mmhij{x,ip-T],v*,aJ[Dxip + r]']) {t,x) > e 



(3.6) 



Let introduce classically {tk,Xk)k a sequence valued in [0,T) x M'^ satisfying {tf^,Xk) —>■ {T,xo) 
and Vi{tk,Xk) — > v*{T,xo). Pick 5 > such that [tk,T] x B{xk,6) C O for fc large enough, and 
introduce the modified test function ip^ given by 

|2 



ip^it^x) := ip{t,x)+[C 



\X X]^ I 



X X]^ 

6 



where <(< T A 6r, (f> is a regular function in C (M ) such that 4>\b{o,i) = 0> 4>\b{o,iY > 0' 
lim|^,|^oo ~ Cfc > is a constant to be determined precisely later on. We deduce 

from ()3.5p that {v* — ip''){t,x) < for {t,x) G [tk,T] x dB{xk,S). Choosing Ck large enough, 
the particular form of the function (p leads to 

c 



{v*-ip'){t,x) < 



for {t,x) £B{xk,dy X [tfc,r]. 



(3.7) 



Thanks to the y/T — t term in the modified test function ip^, we deduce that 
dip'' 



9t 



(t,x) > 0, (3.^ 



for any (t, x, 7]) G [tfc, T) x B{xk, S) x (— T + (, T) and large enough. Choose now ?y < T A | A e 
and introduce the stopping time 6^ := inf {s > ; X|'= ^ B{xk, S) or Jl*^ 7^ } A T , where 
Gk '■= {tk^hXk)- Let finally consider the process (y'^,Z'^', U'^,K'') given on [tfci^fc] by 



^ s 



J s — 

[[v*{s,Xt^)-[^\s,Xt'')-,i 



rs 



dip'' 
~dt 



i6J 



+ Z:^'- + /,e, (., (^* + [ipk - - t;;ejl .^,e,),-ej, ) 



(r,X^'=)(ir 



Lp'^ — rj — Vfk 



One easily checks from ([32])- (1221) -([Ml) that (F'^, Z'', U'',K'') is solution to 



= i0k,x'g^) + / /,e, (x,^s y,, + Ur, Zr)dr - • dw,. 



Ur{j)l-i{dr,dj)+Ke^-Kr 



s Jl 



on [tk,0k], together with the constraint hj^k AX^'' ,Y^- ,Y^- + Ur{j),Zr) > a.e., j G X. Since 
(y^'^, Z'^'^, f/*^*:, K^'') is a minimal solution to this constrained BSDE with jumps and we deduce 

f''{tk,Xk) - V = ip{tk,Xk) + \JT -tk-rj > Vi{tk,Xk), for all A; large enough . 

Letting k go to infinity, this contradicts (j3.5|) and concludes the proof. □ 

Remark 3.2. The main drawback of this representation is the necessity of Assumption (HI). 
Following similar arguments as Proposition 6.3 in [TT], observe that it is satisfied whenever there 
exist a Lipschitz vector function {wi)i^j G [C^(M"')]'^ supersolution to (j3.2p satisfying a linear 
growth condition, and a constant C > such that Owi + /(., {wj)j^j, aj Dwi) < C on M'^, i G X. 
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3.3 A comparison argument 

In this section, we provide sufficient conditions cliaracterizing the value function v as the unique 
viscosity solution to (jl.l|) - (j3.2|) . This gives in particular the continuity of v, leading to the strong 
convergence by penalization and the minimality condition, presented in Section [2j The proof re- 
lies as usual on a comparison argument, which holds under the following additional assumptions. 

(H2) The following holds: 

(i) For any i £ 2, fi is convex in {{yj)j^x, z) and increasing in m. 

(ii) For any i, j € T, /ijj is concave in (yj, yj, z) and decreasing in yj. 

(iii) There exists a nonnegative vector function (Aj)jgj G [C'^{W^)Y' and a positive constant p 
such that, for all i £ I, Ai > gi, \im^,j.^^^ i+\x\ ~ ^^^^ have : 

CAi + fi{., (Aj)jgj, crjD^Ai) < pAi and min hij{.,Ai,Aj,ajD^Ai) > . 

An example where Assumption (H2) holds is given for the case of optimal switching in [2]. 

Remark 3.3. As in Bouchard [2], (iii) allows to construct a nice strict supersolution to (II. ip 
allowing to control solutions to (jl.ip - (j3.2p by convex perturbations. Following the approach of 
[llj . the general form of / and h forces us to add the extra convexity assumptions (i) and (ii). 

Theorem 3.2. Let (H0)-(H1)-(H2) hold. Then, for any U Isc (resp. V use) viscosity super- 
solution (resp. subsolution) to satisfying [\U\ + |l^|](t,2;) < C(l + |2;|) on [0,T] x W^, 
we have Ui > Vi on [0, T] x M.'^, i £ I. In particular, v is continuous and the unique viscosity 
solution to (ll.ip - ()3.2p satisfying a linear growth condition. 

We omit the proof of this comparison theorem which is a natural extension of Theorem 4.1 
in [Tl]. Following the arguments of the proof of Proposition 3.3 in [14], v still interprets as the 
minimal viscosity solution to (|l.ip - (|3.2p in the class of functions with linear growth, whenever 
only a comparison theorem for the IPDE (|3.ip holds. 

4 Numerical issues 

The numerical resolution of systems of variational inequalities of the form (jl.ip - (jl.2p usually relies 
on the use of iterated free boundary. We first solve the system without boundary condition and 
consider recursively the system constrained by the boundary condition coming from the previous 
iteration. In a switching problem, we constrain the solution associated to n + 1 possible switches 
by the obstacle built with the solution where only n switches are allowed. Such a numerical 
approach is computationally demanding. We detail here a natural convergent algorithm based 
on the approximation of the solution to the corresponding constrained BSDE with jump (jl.Sp - 
(jl.6p . We combine a penalization procedure with the discrete time scheme studied by [3] and 
the statistical estimation projection presented in [?]• Thanks to the previous Feynman-Kac 
representation, this gives rise to a convergent probabilistic algorithm solving coupled systems of 
variational inequalities. 
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Let fix an initial condition e £ E and omit it in the notations for ease of presentation. 
Suppose tliat (H0)-(H1)-(H2) holds. The algorithm divides in three steps. 
Step 1. Approximation by penalization. We first approach the constrained BSDE with 
jump (jl.5p - (jl.6|) by its penalized version (j2.ip characterized by a driver /":=/ — as in 

Section [2.21 We deduce from Proposition 12. II that the penalization error converges to as n goes 
to infinity, see (j2.2p . 



Step 2. Time discretization. Observe that the pure Jump process / can be simulated 
perfectly and denote by (t;) its jump times on [0,T]. Let introduce the Euler time scheme 
approximation of the forward process X defined on the concatenation (s/)/ of the regular 
time grid {t^ ■= kh, k = 1, . . . , T/h} with the jumps (r^) of /: 



X^ = Xo and X^^^^ := X^^ + bj^^ {X^,){si+i - si) + aj^^ {X^,)[Ws,^, - Ws,]. 

We deduce an approximation Yj,' of Yf at maturity given by gij.{X!^). The penalized BSDE 
(j2.ip can now be discretized by an extension of the scheme exposed in [3]. An approximation of 
Y^ at time is computed recursively following the backward scheme for k = T/h — 1, • • • , : 



Y, 



k 

•n,h 



^tk+iy^^ik 



E 



tk 



\^n,h p-{{tk,tk+i]x{i}) 
^tk+i A(i) 



i e J 

n,h ryn,h TT'n,h\ 



(4.1) 



ynA rtk+^ jn (j^h M ^n^h 

*fe + l Jtk •'Is^ tfc ' tfc + l' tfe ' tk ■' 



where E^^. denotes the conditional expectation with respect to Gtk ■ Following the arguments of 
Section 2.5 in [3] and identifying (Y"^'^, Z"''*, [/"''^) as a constant by part process on each interval 
{tk,tk+i], we verify the convergence of this time-discretization approximation : 



\Y^ — Y 



n,h I 



+ \\Z'' 



u 



n,h I 



0, 



n G 



(4.2) 



Step 3. Approximation of the conditional expectations. The last step consists in esti- 
mating the conditional expectation operators E^^. arising in ()4.ip . We adopt here the approach 
of Longstaff-Schwarz generalized by [7] relying on least square regressions. 

Fix G N and simulate N independent copies of the Brownian increments (VF/^^^ ~^tk )o<k<T/h 
and the poisson measure {jl^ {{tk,tk+i] x X)(; 



o<k<T/h- For each simulation j < N, define /j and 
^ the trajectories of / and X^. By induction, one can easily verify the Markov property of 
the process (y"=^^ [/",/*) defined in (|iT]) : 



X 



-\rn,h n,h / r v 



ti)^ 



Z^' = a!^\l,^,X^, U^' = l/^\h,,X^, 



for some deterministic functions {a^'^ ^h^'^ ,c^'^)k<n- The idea is to approximate these functions 
using Ordinary Least Square (OLS) estimators. Given L S N, we introduce a collection of basis 
functions {o-f ,hf ,cf)i<i<L of M x M*^ x W^. For each trajectory j < N, define the associated 



terminal value given by Y- 



n,h,L,N 
3,tn 



/vh,N\ TIT 1 j2 -1 / ryn,h,L,N jTn,h,L,N\ 

^ikt ^-^i^tj- Now we define recursively (Z^. '^ ,C/^.^^' ) 
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backward in time for k = T/h — 1, ■ ■ ■ , 0, by computing OLS approximations as follows: 
(/^, ,/:^i)W.-arg^mm^^^2^|^n^^^ 2^ (^.A' J ' 

jr = l / = 1 

for i £ I, leading to the approximation 

L L 
It remains to introduce (71, • • • ,7/,) the minimizer of the mean square error 

j=i -^'fe 1=1 

in order to deduce the OLS approximation Y^l^'^'^ := ^21^=1 ^i^f ^-^ftk'> -^^t^) ■ 

We refer to [7] for the control of the statistical error due to the approximation of the condi- 
tional expectation operators by OLS projections, and, by extension, 

UYn,h_Yn,hL,Nu .U2n,h_^n,h,L,Nu ^ Hljn^h _ ^n,h,L,N u _^ , n G N , h > 0. (4.3) 

The convergence of the algorithm follows from (j2.2p , (j4.2p and (j4.3p . The derivation of a conver- 
gence rate requires precisions on the influence of n on the discretization and statistical errors, as 
well as a control of the penalization error. This challenging point is left for further research. 
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